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Abstract: Based on the calculation of survival probabilities, we suggest a procedure to assess the value 
pf Gsp, the triple Pomeron 'bare' coupling constant, by comparing the large rapidity gap single high mass 
diffraction data in proton-proton scattering and J/^' photo and DIS production. For p-p scattering the 
fcalculation in a three amplitude rescattering eikonal model, predicts the survival probability to be an order 
Q_|Df magnitude smaller than for the two amplitude case. The calculation of the survival probabilities for 
Q_^photo and DIS J/^ production are made in a dedicated model. In this process we show that, even though 
^ its survival probability is considerably larger than in p-p scattering, its value is below unity and cannot be 
• • neglected in the data analysis. We argue that, regardless of the uncertainties in the suggested procedures, 
the outcome is important, both with regards to a realistic estimate of G^p, and the survival probabilities 
^,e,eva„»oLHCe.peri.ne.. 
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1. Introduction 

A large rapidity gap (LRG) process is defined as one where no hadrons are produced in a sufficiently large 
rapidity interval. Diffractive LRG are assumed to be produced by the exchange of a color singlet object 
with quantum numbers of the vacuum, which we will refer to as the Pomeron. We wish to estimate the 
probability that a LRG, which occurs in diffractive events, survives rescattering effects which populate the 
gap with secondary particles coming from the underlying event. 

At high energies, elastic and inelastic diffractive processes account for about 40% of the total p-p (p-p) 
cross section. We would like to remind the reader that: 

1. The small t behavior of the scattering amplitude is determined, mostly, by the large impact parameter 
b values. 

2. The survival probability < jS'p > (denoted S"^) of a diffractive LRG is obtained from a normalized 
convolution of the b-space diffractive amplitude squared and e"^^*'*). 0(s,6) is the optical density, 
also called opacity. Consequently, S'^ decreases with increasing energy due to the growth with energy 
of the interaction input opacity. 
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3. is not only dependent on the probability of the initial state to survive, but is also sensitive to the 
spatial distribution of the partons inside the incoming hadrons and, thus, on the dynamics of the 
whole diffractive part of the scattering matrix. 

4. S*^, at a given energy, is not universal. It depends on the particular diffractive subprocess, as well 
as the kinematic configurations. It also depends on the nature of the color singlet (P, W/Z or 7) 
exchange which is responsible for the LRG. 

Historically, both Dokshitzer et al. [1] and Bjorken [2] suggested utilizing a LRG as a signature for 
Higgs production originating from a W-W fusion sub-process, in hadron-hadron collisions. It turns out that 
LRG processes give a unique opportunity to measure the high energy asymptotic behavior of the amplitudes 
at short distances, where one can calculate the amplitudes using methods developed for perturbative QCD 
(pQCD). Consider a typical LRG process - the production of two jets with large transverse momenta 
Pti ~ —Pt2 ^ ^J'■, with a LRG between the two jets, is a typical mass scale of the soft interactions. 

p{l) + p{2) — > Mi[hadrons + jeti{yi,pti)] (1.1) 
+ LRG[l^y = \yi - y2\] + M2[hadrons + 3et2{y2,Pt2)]- 

yi and y2 are the rapidities of the jets and Ay ^ 1. The production of two hard jets with a LRG is 
initiated by the exchange of a hard Pomeron. We define Fg to be the ratio between the cross section due to 
the above Pomeron exchange, and the inclusive inelastic cross section with the same final state generated 
by gluon exchange. In QGD we do not expect this ratio to decrease as a function of the rapidity gap 
Ay. For a BFKL Pomeron [3], we expect an increase once Ay ^ 1. Using a simple QCD model for the 
Pomeron, in which it is approximated by two gluon exchange [4], Bjorken [2] gave the first estimate for 
Fg ~ 0.15, which is unexpectedly large. 

As noted by Bjorken [2] and GLM [5], one does not measure Fg directly in a LRG experiment. The 
experimentally measured ratio between the number of events with a LRG, and the number of similar events 
without a LRG is not equal to Fg^ but, has to be modified by an extra suppression factor which we call 
the LRG survival probability, 

fgap=<\S\^> >iFg. (1.2) 

The appearance of 5^ in Eq. ( |1.2| ) has a very simple physical interpretation. It is the probability that the 
LRG due to Pomeron exchange, will not be filled by the produced particles (partons and/or hadrons) from 
the rescattering of the spectator partons, or from the emission of bremsstrahlung gluons coming from the 
partons, or the hard Pomeron, taking part in the hard interaction. 

<| S P> = <| Sbremsstrahlungi^y = | ^1 "2/2 1) P> X <| S gpectator { S ) |^> , (1.3) 

where s denotes the total cm. energy squared. 

• ^bremsstrahiung(^y) ^^^^ Calculated in pQCD [6]. It depends on the kinematics of each specific 
process, and on the value of the LRG . This suppression is commonly included in the calculation of 
the hard LRG sub process. 
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• To calculate S^p^^^^^^jX^) '^^ need to find the probability that all partons with rapidity yi > yi in 
the first hadron, and all partons with yj < y2 in the second hadron, do not interact inelastically 
and, hence, do not produce additional hadrons within the LRG interval. This is a difficult problem, 
since not only partons at short distances contribute to such a calculation, but also partons at long 
distances, for which the pQCD approach is not valid. Many attempts have been made to estimate 
'^spectator [2,5, 7^11], but a Unique solution, has still not been found. 

An obvious check of the above is to compare the calculated values of Sgp^^i^fgj.{s) obtained in different 
models for different reactions. The Durham group [12] recently suggested a very interesting procedure, 
proposing to extract G^p, the triple Pomeron vertex coupling, utilizing the measurement of large mass 
diffraction dissociation in the reaction 

7*(Q^^Bj) + P ^ J/^ + [LRG] + X{M^ > rup). (1.4) 

The cross section of this process can be described by the Mueller diagram (Fig. |l]). It is initiated by the 
charm component of the photon which has a small absorptive cross section, since its interaction stems from 
short distances (r oc l/rric, where rric is the mass of the charm quark). Thus, the probability for additional 
rescatterings (Fig. |2|) is relatively small, resulting in a high survival probability. This is to be compared 
with the corresponding high mass diffraction in an hadronic p-p (p-p) reaction (Fig. |^, for which we expect 
the rescatterings (Fig. ^ to be significant, resulting in a small survival probability [13]. It is, therefore, very 
probable that present extractions [14, 15] of the 3P coupling are underestimated. Comparing the values of 
Gsp obtained in the above two channels, taking into account their (different) survival probabilities, leads 
to a more reliable measure of the 3P coupling, and provides a check of the various theoretical estimates of 
the survival probabilities. 




Figure 1: The general diagram for difFractive pro- 
duction of large masses in "f*-p collisions at high en- 
ergy. Pomcrons are denoted by the zigzag lines. The 
vertical lower Pomeron is soft. The two Pomerons 
coupled to the photon vertices are hard. 



Figure 2: The general diagram for calculating 
the survival probability for difFractive production of 
large masses in j*-p collisions at high energy. The 
Pomerons identification is identical to the previous 
figure. 



The purpose of this paper is to re-examine the procedure to extract the value of Gsp, the triple 
Pomeron 'bare' coupling constant. As a by-product we assess the stability of the survival probabilities 
obtained in two and three amplitude eikonal rescattering models. The following topics are addressed: 
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Figure 3: The general diagram for diffractive pro- 
duction of large masses in p-p collisions at high en- 
ergy. The zigzag lines denote soft Pomerons. 



Figure 4: The general diagram for calculating the 
survival probability for diffractive production of large 
masses in p-p collisions at high energy. The zigzag 
lines denote soft Pomerons. 



1. We estimate the survival probability values for the triple Pomeron vertex in p-p high mass single 
diffractive (SD) collisions. 

2. We investigate the difference between the survival probabilities associated with high mass SD and 
the inclusive SD channels. 

3. The above calculations are carried out in the two channel GLM model [16-18] in which we compare 
the output of its two and three amplitude representations, so as to assess the reliability of our results. 



4. We calculate the survival probability for the reaction given by Eq. (|l.4| ). Note that in an hadronic 
high mass SD reaction the 3P coupling consists of three soft Pomerons. Eq. (|l.4| ) is a hard process 
in which two hard Pomerons couple to a soft one. It is possible, but not necessary, that the above 
3P couplings are equal. 

The plan of our paper is as follows: In the next section we calculate Sgp^^fatori^) foi' high mass 
diffraction in p-p scattering in two [16] and three [17] amplitude models for the soft interactions. The 
models are specified and their outputs compared. In section 3 we present our estimates for the survival 



probability of the process of Eq. (L4). These calculations, carried out in a dedicated model [19], show that 
even though the J/^ survival probabilities are significantly larger than those calculated for p-p, they cannot 
be neglected. We note that, the reliability and accuracy of the J/^ calculations are considerably better 
than in the p-p channel. In the conclusions, we summarize our main results and specify some remaining 
problems. 



2. Survival probability for the triple Pomeron vertex in proton-proton collisions 
2.1 Survival probability in the eikonal model 

The cross section for diffractive dissociation in the region of large M can be viewed as a Mueller diagram 
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(Fig. 1) which can be rewritten in terms of the triple Pomeron vertex (see Ref. [14]). We denote this cross 
section a^^ and its corresponding survival probability, at a given M^, S^p{M'^). 



where g{t) describe the vertex of Pomeron-proton interaction, and G^p stands for the triple Pomeron 
vertex. However, this diagram does not take into account the possibility of additional rescatterings of the 
interacting particles shown in Fig. ^. The result can be written as 

^ dtdM^ " '^3p(^n (^^j ■ (2.2) 

The survival probability factor S'|p(M^) is defined ^ as 

r kt M2 -^^^ {Fig. H) 

^-<-^' ^ k^i^is^ ■ 

The easiest way to calculate the diagram of Fig. ^ is to first transform the diagram of Fig. ^ to impact 
parameter space. This is done by introducing the momentum q along the lowest Pomeron in Fig. |3|. In this 

case, 

T{s,M'-g) ^ J d'hM'^^^{F^g.^ (2.4) 

2 g,ik^) g,{{k - q)l) g^jq^) G^p{kl {k - q)lq^) / s x ap(fc?) + «p((fc-,-0?) - 2 / \ 

16^^ \W) 

Prom Eq. (^3) we find the form of this amplitude in impact parameter space to be 

T{s,M^-h) ^ I -0^ A{s,M^;q). (2.5) 

Using a linear approximation for the Pomeron trajectory and a Gaussian form for all vertices 

apit) = l + Aa'pt, gp{k^)=gp{d)e-^^^\ G3p(A:i, ^2, A^g) = G3p(0,0,0) e'^-^'^i + '^' + ^i) , (2.6) 
we obtain 

rr ( A#2 ,^ 53P ic\ I \ ^ ( d{y) [d{0 + d{0] /o 7^ 

S = ln(A/Vso), J = ln(s/A/'), 93P s O3p(0,0,0)/9p(0), ,,(■;,) = -SM^ e^f, (2.8) 



S^p{M ) denotes the high mass SD survival probabihty, at a given M , and is identical to Sgp^ctatoris) for this specific 
SD reaction. 
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where ^ 

R\y) = 2Rl + 2rl + Aa'py and d{y) = (2.9) 

Tg denotes the radius of the triple Pomeron vertex. In the following we take [20] Tq = 0.5 GeF"^. 

Using Eq. (2.7), the expression for the survival probability (see Eq. ( |2.3| )) in a simple eikonal model, 
accounting for the rescattering corrections, can be written as 

„ , „ r (i26r(s,M2;6) exp(-fi(C + y;6)) , , ^2 

Slp{M^) = ^ fd^Tis M^;b) ' n{C + y;b) = i^ppe (2.10) 

""PP = A ^^^^""''^ Rpp{^ + y) = ARIp + Aa'pi^y). (2.11) 

■KRpp[i + y) 

2.2 Two channel models: main ideas and formulae 

In the eikonal model only elastic rescatterings have been taken into account. Two channel eikonal models 
have been developed so as to also include rescatterings through diffractive dissociation (see Refs. [16-18] 
and references therein). In this formalism, diffractively produced hadrons at a given vertex are considered 
as a single hadronic state described by the wave function which is orthonormal to the wave function 
of the hadron (proton in the case of interest), < o > = 0. 

Introducing two wave functions that diagonalize the 2x2 interaction matrix T 

Ai^k = < *i ^fc|T|^i/ -^k' > = Ai^k 5k,k', (2.12) 

In our past publications we referred to the GLM eikonal models according to the number of the rescattering 
channels considered, i.e. elastic [20], elastic + SD [16] and elastic + SD + DD [17]. In retrospect, we 
consider it more appropriate to define these models according to the dimensionality of their base. We, 
therefore, call the above a two channel model, making the distinction between its two and three amplitude 
representations. 

We can rewrite the amplitude A^^k in a form that satisfies the unitarity constraints 

AAs.b) = i fl - exp (_^hdhR)) . (2.13) 



In this formalism we have 

Gf^.is, b) = 1- exp (-0,,fc(s, b)) . (2.14) 

G*" is the probability for all inelastic interactions in the scattering of particle i off particle k. From 
Eq. ( 2.14| ) we deduce that the probability that the initial projectiles reach the interaction unchanged, 
regardless of the initial state rescatterings, is exp (—0,^^(5, 6)). 

In this representation the observed states can be written in the form 

= a^^i + P^2, = + , where + = 1. (2.15) 
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The obvious generalization of Eq. (|2.7D is 

T (s, M^; 6) = ^ < p\l >2 < p|A; > T^'^ (s, 6) < p|A; >< p\i (2.16) 

where < p\l > = a and < p|2 > = /3. 

{',M = ^ My) ,,(^) ^ ^ exp ^^^^^ ^ ^^^^^ >, j , (2.17) 

where S'sp = G3p/gi{0) and 

My) = ^-M^Me^y and = 2RI, + 2rl + 4a'py, 4(2/) = (2.18) 
The numerator of Eq. ( p.lO| ) is written in this representation as 



(fb <p\l>'^<p\k> e — Tfj' (s, M^; b) e — < p\k >< p\i 
For J7j^fc(s, 6) we take 



^^,k{^,b) = u,,ke , (2.20) 

where 

= ^^^f^f-)^, and RUC + y) = 2Rl + 2Rl, + Aa'AC + y). (2.21) 

In the following we denote g'^^ = gi{0)gk{0)- 

In our calculations we take [16] i?g ^ as a free parameter while i?Q 2 = 0. The survival probability can 
be calculated as the ratio 

S2 .^2^ _ IdHN{C,y;b) 



where 



and 



+ a^/3^r2^'^(6)e-^^>^W + 2 a' T^'\b) e-^^^^^^^^^ + /^^ T^^'^^fe) e-^^.^('') , (2.23) 



D(e,y;6) = a''Tl'\b) + 2a^fi^Tl^\b) + a^p^T^'^b) 

+ a^P^T^'\b) + 2a2/34r2^'2(5) + T^^\b) . (2.24) 
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For completeness we also present the integrated cross sections of the diffractive channels in the two 
channel model, together with the corresponding elastic and total cross sections. The amplitudes for the 
elastic and the diffractive channels have the following form [16, 17] 

a^i {s,h) = i il - e 2 — - 2a^ [3^ e 2 — - (3^ e 2 — J , (2.25) 

(^1 i(s,b) fil 2(s,6) f22 2 

a^e ^ — - (a^ - I3^)e ^ — - /^^ e ^ — ) , (2.26) 

add{s,b) = ia^ (3^[ -e 2 — + 2e 2 — - e 2 — . (2.27) 



Using Eq. (|2.20|) and Eq. ( 2.21 ), as well as the general cross section formulae we get 



atot{s) = 2 J <fbaeiis,b), aeiis) = J cf b \aei (s , b) {"^ , (2.28) 

cTsdis) = J dH\asdis,b)\^, addis) = j b\aM{s,b)? ■ (2-29) 

It is instructive to present the calculation for the diffractive channels in the form of a survival probabil- 
ity, which we define as the ratio of the output corrected diffractive cross section to the input non corrected 
cross section. 

2 _ / d?b\asd{sM? o^^ 
~ /dn|ay.,6)|2' ^'-"""^ 

where 

«L {s, b) = '-^{- a' ni^iis, b) + [a^ - 0") ^1,2(5, b) + /J^ (^2,2(5, b)) . (2.31) 
We will discuss the results and interpretation of our calculations in the next subsection. 

2.3 Three models used for fitting the experimental data 

To calculate the survival probabilities one needs to specify the opacities VLi^k{s-,b). These are determined 
from a global fit of the experimental soft scattering data. We have used three models based on the general 
formulae given in Eq. ( p. 17 ) - Eq. ( p. 21 ), but with different input assumptions. Note that the above 



global fit has in addition to the Pomeron contribution, also a secondary Regge sector (see Ref. [16]). 
This is necessary as the data base contains many experimental points from lower ISR energies. A study 
of the Pomeron component alone, without a Regge contribution, is not possible at this time, since the 
corresponding high energy sector of the data base is too small to constrain the fitted parameters. The 
Regge parameters are not quoted in this paper and will be discussed in detail in a forthcoming publication. 



- 8 - 



2.3.1 Two amplitude model 



In this model, denoted by Model A [16] we assume that the double diffraction cross section is negligible, 
and we take a^d in Eq- (|2.27l ) to be zero. This allows us to express ^2,2 in terms of Jli^i and ^^1,2, leading 
to the following formulae (see Refs. [16-18]): 

„.Js.,)^i(,-..J-^)-2l,'..J-^) A-expf-^))V (2,32) 



. ,x ^ / An(s,b)\ f / An(s,b)\ , 

asd{s,b) = -ia(3 exp ^ 1 - exp ^ . (2.33) 



The above is a two amplitude model with two opacities l^i^i and A^l = 0,1^2 — ^1,1 • Following Ref. [16], 
we assume both Qi^i and AO to be Gaussian in b. 




?2 



Note that in this two amplitude model Rj^ is the radius of Ai7(s, b). As we shall see, in the three amplitude 
model Ri 2 is the radius of Oi,2- The radii Rf i and Rf 2 are specified in Eq. (2.21). We have also studied 



a two amplitude model in which both l^i^i and J7i^2 are Gaussian in b. The output obtained in this two 
amplitude model is compatible with the output of Model A. This is a consequence of the fit, discussed 
below, in which consistently Oi^2 ^ 

2.3.2 Three amplitudes models 

In the three amplitude models we do not make any assumptions regarding the value of the double diffraction 
cross sections [21] which are contained in our data base. We use Eq. ( 2.2C| ) and Eq. ( |2.21| ) to parameterize 



the three independent opacities: ^^1,1, ^^1,2 and 1^2,2) which are all taken to be Gaussian in b. For details 
see Ref. [17]. Assuming Regge factorization, ^1^2 is determined by ili^i and 0,2^2, see Eq. ( 2.21| )). We 



denote this model B(l). As we shall see. Model B(l) does not reproduce the data well. We have, thus, 



also examined a non factorizable model, denoted B(2), in which ui^2 of Eq. ( 2.21 ) has been replaced by the 
expression 

-1. = , (-)\ (2.36) 

where 5^ 2 is a free parameter. This additional degree of freedom violates Regge factorization for the input 
Pomeron, but it allows us to describe the experimental data on the double diffraction cross section which 
we failed to fit in model B(l). 
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2.4 Results 



The parameters of Model A , quoted from Ref. [16] are based on a fit to 55 experimental data points 
base which include the p-p and p-p total cross sections, integrated elastic cross sections, integrated single 
diffraction cross sections, and the forward slope of the elastic cross section in the ISR-Tevatron energy 
range. As stated, we neglected the (very few) reported DD cross sections. The fitted parameters of Model 
A are listed in Table 1 with a corresponding / {d.o.f) of 1.50. 

The fit to Models B(l) and B(2) are based on the formulae given in the previous subsection and in 
Refs. [17, 18]. The new data base includes the data used for Model A, plus 5 double diffraction cross 
sections points [21]. The values of the parameters and associated with B models are also given in 
Table 1. The best fit obtained from Model B(l) has a / {d.o.f) = 2.3, which is unsatisfactory. This is 
manifested not only by the high x^, but also by inability of the model to reproduce the Udd points, which 
are severely underestimated. We thus consider Model B(2) in which the couplings of the three amplitudes 
are independent. This model produces a much better fit with / {d.o.f) = 1.25, as well as reproducing 
the experimental results for add- The cross section predictions of Model B(2) are shown in Fig. ^ The 
parameters associated with Models A, B(l) and B(2) were used to calculate the survival probability for 
two different soft SD final states. S'^p corresponds to single diffraction dissociation in the high mass region 
(see Eq. ( |2.22 )). We note that S'|p is weakly dependent on large - mass of produced hadrons. is the 
survival probability corresponding to the entire region of the produced diffractive mass (see Eq. ( |23oD ). 
The calculated survival probabilities are presented in Fig. ^. We note a significant difference between the 
various outputs which will be discussed below. 

The following are a few qualitative characteristic remarks and comments: 



1. Clearly, our analysis favors the non factorizable input of Models A and B(2) over the factorizable 
input of Model B(l), which is theoretically more appealing. Technically, the factorization breaking 
is induced by the output results, in which §2 is considerably larger than gi resulting in Ai^2 S> 

As we have shown in Ref. [17], this is a consequence of the striking experimental observation that 
in p-p (and p-p) scattering Rd = {uei + (Jgd + crdd)/(^tot ~ 0.4. This is different from Rei = Oei/cJtot 
which rises monotonically with energy. The analysis made in Ref. [17] found that to be consistent 
with the experimental behaviour of Rd and R^i requires that g2 > 10 (71. The best value obtained is 
g2 = \/300gi which is very close to the values found in Model B(2). 

2. We observe a general systematic behaviour in which the 3P and SD survival probabilities become 
smaller when we add an amplitude to the initial state rescattering chain. In the context of this paper 
we find that S'^ (Model A) > (Model B) for both 3P and SD channels. 

We trace this dramatic difference to the relaxation of the constraints on the initial rescatterings. 
Allowing an additional amplitude, enables increased screening of the input amplitude. 

3. An important observation, correlated to the above, is that the unscreened input cross sections of 
Model B(2) are considerably larger than the unscreened input cross sections of Model A. This is a 



- 10 - 



■ Total, elastic and diffractive cross sections 




10^ 10"* W(GeV) 10^ 10 



Figure 5: Total, elastic and diffractive dissociation cross sections in model B(2). 



Model 


A 


P 


^0;1,1 


/ 

Qp 


all 


92,2 


gi 


91,2 


A 


0.126 


0.464 


16.34 GeV-^ 


0.200 GeV-^ 


12.99 GeV~'^ 


N/A 


145.6 GeV'^ 


N/A 


B(l) 


0.150 


0.526 


20.80 


0.184 GeV-^ 


4.84 GeV-^ 


4006.9 GeV-^ 


N/A 


139.3 GeV-'^ 


B(2) 


0.150 


0.776 


20.83 GeV'^ 


0.173 GeV'^ 


9.22 GeV-^ 


3503.5 GeV-^ 


N/A 


6.5 GeV-^ 



Table 1: Fitted parameters for Models A, B(l) and B(2). Rlj^ ^ ^o-a = 0-5^o;i,i' ^0;2.2 = 0- 



consequence of the large difference between tlie fitted values of g'l 2 of Model B(2) and g"^ of Model 
A, and, also, the different values of A (the Pomeron intercept) of these models. Since these models 
provide compatible good reproductions of the fitted data base, we conclude that the large difference 
in the unscreened cross sections of Models A and B(2) are compensated by the reciprocal difference 
in the corresponding survival probabilities. 

4. The survival probabilities calculated in our two amplitude Model A are in agreement with those 
calculated by Khoze, Martin and Ryskin in their two channel model [12], which differs from ours. 
Our present observation that the three amplitude Model B results in considerably smaller S'^ values, 
implies that the presumed consistency between most of the published survival probability outputs 
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Figure 6: Survival probabilities for single diffraction p-p collision. S'^p (solid line) denotes the survival probability 
for high mass diffraction dissociation (see Eq. ( ^.22 )). S^^ (dotted line) is the survival probability for diffractive 



dissociation in the entire kinematic region (see Eq. (|2.30| ). The upper curves (solid and dotted) refer to Model A [16]. 
The lower curves relate to Model B of Ref. [17]. 

(for details see Ref. [18]) should be carefully re-examined using more robust models. 

5. 5|p are consistently higher than 5^^. Note that the input of crsd(high mass), described by the triple 
Pomeron diagram, behaves as s^. On the other hand o"srf(low mass), for which the triple Pomeron 
diagram is not applicable, behaves as s^^. 

Diffractive processes are very important at cosmic ray energies. Kama et al. [22] have shown that 
diffractive p-p interactions play a crucial role in understanding the spectrum of galactic gamma-rays that 
come predominantly from vr'^ 77. The inclusion of diffractive processes makes the gamma ray spectrum 
harder, and when this is included together with the assumption of Feynman scaling violations, one can 
explain about half of the " GeV Excess" . The magnitude of the diffractive cross sections is a crucial element 
of this and similar studies. A review of our results at cosmic ray energies will be published elsewhere. 

3. Survival probability for triple Pomeron vertex in photo and DIS production 



In this section we calculate the survival probabilities for high mass diffraction in the reaction of Eq. (1.4). 



From Fig. ^ one can see that we need the following ingredients to make these estimates: i) the amplitude 
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for the interaction of a colorless dipole with the target.; and ii) the description of J/^' production with no 
initial state interactions with the target, shown in Fig. |l]. 

For the scattering dipole amplitude we take a model developed by one of us [19]. This model is based 
on the solution of a generating functional [23-25], with an additional assumption that the dipoles do not 
change their sizes during the interaction. The amplitude is equal to 



N{Y = ln(l/x);r, b) 



niY;r,b) 
1 + niY;r,b) ' 



where 



niY;r,b) = ^j^xG{xo,f^')xG{^,i,' = S + /"o ) Sib). 



S{b) is the proton b-profile. 



7ri?2 R 



R 



(3.1) 



(3.2) 



(3.3) 



The saturation scale Q'^{Y) is defined from the condition 



n{Y;rs,b) = 1 with = i/Ql 



(3.4) 



The gluon structure function xG{x,fj?) satisfies the DGLAP evolution equation with the initial con- 
dition xG{x,Qq) = A/x'^" at Qq = 1 GeV^. In other words, Eq. p.2| ) describes the contribution of the 



hard Pomeron that can be calculated in pQCD. All parameters in Eq. ( |3.2D have been found by fitting to 
the data on F2. The fit is good and has x^/d-o.f. = 1.07 [19]. The values of the fitted parameters, are 
compatible with the parameter values obtained in other competing models [26]. Eq. ( |3.3D is the Fourier 
transform of the electro-magnetic form factor of the proton in t-space. Eq. (|3.1| ) is quite different from the 
eikonal approximation that has been used in other models, and has a form which is typical for the 'fan' 
diagrams which are summed in the mean field approximation (MFA). 

One can see directly from Eq. (|3.1|) and Eq. (p^) that 



N{Y = ln(l/x);r, b) 



r — ► 00; X fixed 



N{Y = ln(l/x);r, b) 



■- 0; r fixed 



1. 



(3.5) 



In 



Fig. ^ we show how Eq. (3.1) fits the experimental data. 

re can write the formula for Fig. 

J/^ + M{Fig. |lD) = GsP (m2)^^ I dp b exp (^--^^ x (3.6) 

dz(f r'^^*{r,z,Q'^)A{r,xp)^jl^{r,z) x j dz' (f r'^^*{r' z' ,Q'^)A{r\xp)^ ji^{r' ,z' 



Using Eq. ( |3.2D we can 

CT (7* + p 
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= 0.25 = 0.5 = 2 




Figure 7: Examples of the fit with Eq. ( |3.l|) . 



where 

A{r,xp;b) = j cfbn{Y;r,b). (3.7) 

In Eq. ( p.6p we assume that i) that the hard Pomeron in the upper legs of Fig. ^ does not depend on the 
impact parameter, ii) The triple Pomeron vertex (HP — P — HP) couphng is the same as the triple soft 
Pomeron vertex {P — P — P) coupling. HP denotes the hard Pomeron, while P stands for the soft Pomeron. 
R is the radius of the soft interaction and it was taken to be equal to B? = (12 + In(MVso)) GeF"^ 
with So = IGeV"^. 

The product of the wave functions is taken as [27] 

|£ (^ci^,(c)-^i^2(C))+"^c (^K2(C)-Ci^i(C))}, (3.8) 



^,/,(.,. = -) X ^,,,(.;Q^) = 



^f^i^2(C)-Ci^i(C))^ (3.9) 



2 V 2 
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C = ar, Ki (where i=l,2) are the modified Bessel functions, Fee = 5.26 KeV is the leptonic width of J/^'. 

Q2 + m2 



a = z {1 — z) Q + m^, and xp 



(3.10) 



The exponential factor in Eq. ( p.§| ) and Eq. (3.£) which describes the wave function of J/^* meson with 
velocity v, has been discussed in Ref. [28] and references therein. 

For the diagram of Fig. Q we can write the expression that takes into account the possible rescatterings 
before the interaction that produces a J/^' meson, 



a (7* + p ^ J/^ + M{Fig. D) = GsP (M^)^^ J dH exp , ^ , 
dzd'^r^y*{r,z,Q'^) A{r,xp) (1 - N{Y = ln{l/x);r,b)) ^j/^{r,z) 
dz' d^r''^-y*{r,' z',Q^) A{r',xp) (l - N{Y = ln{l/x);r' ,b)) ^'j/vi,(r', z') 



(3.11) 



The survival probability is the ratio of these two equations (Eq. (|3.6|) and Eq. ( 3.11| )) 

a{j*+p^ J/^ + MjFig. |; Eq. (pl|))) 



^3P('//^) 



(7 (7* + p ^ J/l- + M(Fig. |; ^g. (^) ) ' 



(3.12) 



The results of our calculations using Eq. ( 3.121 ) are plotted in Fig. ^. One can see that S^p is smaller 
than 1, not compatible with the assumption of Ref. [12]. This should be taken into account when attempting 
to extract the value of the triple Pomeron vertex from the measurement of the cross section of reaction of 
Eq. Q. 

Taking the survival probabilities corrections into account when extracting the value of G^p, it is 
interesting to study the ratio R = Rp/Rj/^, where 



Rn 



R 



J/* 



a{p + p ^ p + p) 
a{Y+P^ J/^ + M(M > rup)) 



(3.13) 



(3.14) 



Gsp ir+p^ J/^ + p) 

so as to validate the theoretical estimates, as well as checking the sensitivity of G^p to the hardness of the 
coupled Pomerons. 



4. Conclusions 

In this paper we confirm the wide spread expectation that the survival probability for the triple Pomeron 
vertex is very small [12, 13, 15]. However, whereas we find the results of our two amplitude Model A almost 
identical to the results obtained in Ref. [12] (which is also a two amplitude model), our three amplitude 
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Figure 8: The survival probability for the triple Pomeron vertex in the reaction of Eq. ( |l.4| ). 

model estimates are an order of magnitude smaller. This may also influence the 5^ calculations for other 
channels. In particular, it was noted [18] that the various two amplitude calculations of S'^j-,, relevant 
to exclusive central Higgs production, are in remarkable agreement. This optimistic evaluation should be 
carefully re-examined. 

We stress the importance of J /'^ SD photoproduction (Eq. ( |1.4| )). Comparing its high mass diffraction 
data with the corresponding p-p scattering, and correcting both channels with their corresponding survival 
probabilities, we hope to evaluate both the value of G^p and its dependence on the Pomeron's hardness. 

In this context, we emphasize that even though S'|p obtained for Eq. is considerably higher than 
S'^p obtained for p-p scattering, its value is less than unity and should not be neglected. 

The details of our three amplitude model B, including calculations for hard diffraction, and other LRG 
configurations will be published in a forthcoming paper. 
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